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Experimental conclusions from air shower observations on cosmic-ray photons above 
10^® eV are based on the comparison to detailed shower simulations. For the calculations, 
the photonuclear cross-section needs to be extrapolated over several orders of magnitude in 
energy. The uncertainty from the cross-section extrapolation translates into an uncertainty 
of the predicted shower features for primary photons and, thus, into uncertainties for a 
possible data interpretation. After briefly reviewing the current status of ultra-high energy 
photon studies, the impact of the uncertainty of the photonuclear cross-section for shower 
calculations is investigated. Estimates for the uncertainties in the main shower observables 
are provided. Photon discrimination is shown to be possible even for rapidly rising cross- 
sections. When photon-initiated showers are identified, it is argued that the sensitivity of 
photon shower observables to the photonuclear cross-section can in turn be exploited to 
constrain the cross-section at energies not accessible at colliders. 

PACS: 96.40.Pq,96.40.-z,13.85.-t,13.85.Tp 

Key words: Cosmic rays, air shower, photon, cross-section 


1 Introduction 

Photons around and above 10^® eV might provide a key to understanding the 
origin of cosmic rays. Substantial fluxes of these ultra-high energy (UHE) photons 
are predicted in top-down models of cosmic-ray origin. In addition, UHE photons 
are produced by the GZK process of resonant photoproduction of pions m , in 
analogy to GZK neutrinos. 

Experimentally, UHE photons can be discerned from nuclear primaries due to 
differences in the expected shower signatures. So far, no claim of a photon detec¬ 
tion exists and upper limits to UHE photons were set. Any conclusions about UHE 
photons rely on the comparison of data to detailed simulations of photon-induced 
showers. Although these photon showers are dominated by electromagnetic inter¬ 
actions, a source of uncertainty is the photonuclear cross-section, which has to be 
extrapolated over several orders of magnitude in energy from laboratory data for 
calculating showers induced by UHE photons. 

In this work, after giving an overview of the current status of UHE photon 
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Fig. 1. Production of photons by the GZK process; Bottom-up sources were fitted to 
the HiRes spectrum (left) and the AGASA spectrum (right). Plotted is the resulting 
photon fraction, integrated above the primary energy, as a function of the primary energy 
threshold. The ranges shown for the GZK photons come from varying the source and 
background parameters. Figures are taken from [7|. 


studies, both the impact of the uncertainty from the photonuclear cross-section to 
photon shower features and possible prospects for constraining the cross-section 
are investigated. The plan of the paper is as follows. In Section [3 scenarios for 
producing cosmic-ray photons above 10^® eV = 10 EeV are briefly described. Fea¬ 
tures of showers induced by photons are discussed in Section |21 In Sectional the 
experimental situation concerning upper limits to photons is summarized. The role 
of the photonuclear cross-section is investigated in Sectional 

2 Cosmic-ray photons 

Motivated in particular by reports from the AGASA experiment about a pos¬ 
sible continuation of the cosmic-ray energy spectrum without a GZK cutoff, many 
top-down scenarios were proposed P]- In these models, cosmic rays are produced by 
the decay or annihilation of superheavy dark matter or topological defects ^ . A 
common feature of these models is the prediction of a significant flux of photons ar¬ 
riving at Earth. Also in the Z-burst model [S], large photon fractions are predicted. 
Stringent upper limits to photons would disfavour many realizations of these top- 
down models, although some uncertainty exists in calculating the propagation of 
UHE photons [HI- 

Even in “conventional” bottom-up models, in which nuclear particles are accel¬ 
erated in astrophysical sites to highest energies, UHE photons are expected. They 
arise from the GZK process during propagation. Such GZK photons were stud¬ 
ied in detail recently in |7]. In Figure ^ results of these calculations are shown: 
For various assumptions on source parameters (distribution, emission features) and 
on radio backgrounds and magnetic fields important for photon propagation, fits 
were performed both to the HiRes (cutoff-like) [S| and AGASA (no-cutoff-like) |H1 
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cosmic-ray spectra. 

In case of the HiRes spectrum, the range of possible photon fractions extends 
from almost negligible (<0.01%) to values of 1% above 10 EeV and 5% above 
100 EeV. For the completed Pierre Auger Observatory (Northern and Southern 
site) [in] . the latter values would imply event rates of up to ~60 photons/year above 
10 EeV and ~2-3 photons/year above 100 EeV. Such a signal seems detectable; in 
turn, an absence of these GZK photons, and correspondingly an upper limit to 
photons, might translate into constraints for bottom-up scenarios. 

In case of the AGASA spectrum, it turned out that even bottom-up scenarios 
might allow a fit to a no-cutoff spectrum. The GZK photon fraction is then very 
large, however. (The term “GZK photons” refers to the way these photons are 
produced; it does not imply that the observed energy spectrum necessarily exhibits 
a GZK cutoff.) Thus, a connection between the shape of the energy spectrum and 
the photon fraction predicted in bottom-up models exists [ 7 ]. It seems worthwhile 
to note that in turn, a small-enough photon limit (e.g. below 2-3% above 10 EeV or 
below 10-20% above 100 EeV) would then disfavour a no-cutoff-like spectrum for 
the scenarios regarded in | 7 ], unless specific new physics such as Lorentz invariance 
violation is invoked (of course, top-down models would be constrained by the same 
photon limit). 


3 Air showers induced by photons 

In FigureEl the average depth of shower maximum Amax is shown vs. primary 
energy as calculated for different primary particles. The shower maxima of primary 
photons are separated from those of nuclear primaries by ~200 g cm“^ or more 
at 10 EeV. The elongation rate (slope of the curve) is larger for photons even at 
relatively small energies. This elongation rate for photons can be understood within 
the Heitler toy model of shower development m The slope is increased above a 
few EeV with the LPM effect ^21 becoming increasingly important. At even higher 
energy (above 50 EeV, depending on geomagnetic field), UHE photons may convert 
in the geomagnetic field and form a preshower na before entering the atmosphere. 

The directional dependence of the preshower effect is illustrated in Figure |21 
(left). In this example, photons of 100 EeV primary energy do not convert when 
entering almost parallel to the local geomagnetic field lines. For a more perpendic¬ 
ular incidence, in most cases a preshower is formed m 

If a preshower is formed, the energy is distributed among secondary electrons, 
positrons and, mostly, photons. As a specific example, in Figure|21 (right) the energy 
spectrum of preshower particles on top of the atmosphere is shown. As the original 
particle, a primary photon was assumed for the geometry of the 320 EeV Fly’s 
Eye event m On average, about 1400 lower-energy particles are produced. Most 
energy is carried by particles of ~10 EeV energy, i.e. about 1.5 orders of magnitude 
below the original primary energy ini 

The impact of the LPM and preshower effect is shown for the example of the 
highest energy Fly’s Eye event in Figure 01 Assuming primary photons, the longi- 
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Fig. 2. Average depth of shower maximum Xmax versus primary energy simulated for 
primary photons, protons and iron nuclei. For nuclear primaries, calculations for different 
hadronic interaction models are shown. Also shown are experimental data. For references 
to the experiments and interaction models see m 


tudinal profiles are shown for calculations with both effects neglected, with only 
LPM switched on, and for the complete treatment. Due to the extreme energy of 
the Fly’s Eye event, differences between the resulting profiles are enormous. While 
the LPM effect delays the development and increases fluctuations, the preshower 
effect counteracts. 


4 Data 

For the 320 EeV Fly’s Eye event, the observed profile m differs from the pri¬ 
mary photon prediction by about 1.5 standard deviations [mini. The comparison 
is shown in Figure 01 (right). A photon origin of the Fly’s Eye event can not be 
excluded; profiles calculated for nuclear primaries fit the data better, however m 
Limits to photons were set by different experiments. Comparing rates of nearly 
vertical showers to inclined ones, upper limits to the photon fraction of 48% above 
10 EeV and 50% above 40 EeV (95% CL) were deduced from Haverah Park 
data m Based on an analysis of muons in air showers observed by the Akeno 
Giant Air Shower Array (AG AS A), upper limits were estimated to be 28% above 
10 EeV and 67% above 32 EeV (95% CL) [201. In a dedicated analysis of the 
highest-energy AGASA events, an upper limit of 67% above 125 EeV (95% CL) 
was set m These limits came from ground arrays. Recently, using Xmax observed 
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Fig. 3. Preshower features; (Left) Conversion probability versus primary energy for two 
different arrival directions: Small angle (“weak B”) and almost perpendicular (“strong 
B”) to the local magnetic field. Calculation performed for the Auger Southern site ca 
(Right) Spectrum of preshower particles (with and without energy weighting) on top of the 
atmosphere calculated for photons for the geometry of the 320 EeV Fly’s Eye event 




Fig. 4. (Left) Photon shower profiles for the geometry of the 320 EeV Fly’s Eye event 
for the full simulation (solid), when switching off the preshower effect (dotted-dashed), 
and when switching off also the LPM effect (dotted) ca (Right) Observed profiles of the 
320 EeV Fly’s Eye event Util compared to the full simulation for primary photons 1171 . 
Data points are correlated with respect to shifts in atmospheric depth X. 


by fluorescence telescopes in hybrid events, a limit of 26% above 10 EeV (95% 
CL) was obtained by the Pierre Auger Collaboration m These upper limits are 
compared to some model predictions in Figure El 

With increased event statistics, much stronger constraints on the photon flux or 
the discovery of UHE photons can be expected. There exists a minimum possible 
value for an upper limit that ideally could be reached for a given event statistics. 
This is due to the fact that, assuming a fraction of photons in the primary flux, 
a set of rim primaries picked at random is expected to ab initio contain no primary 
photon with probability (1 — Ely)"™. The relation between the minimum possible 
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Fig. 5. Upper limits (95% CL) to cosmic-ray photon fraction derived from Haverah 
Park (HP) 1191 and AGASA (Al) 1201 . (A2) 1211 data and by the Auger Observatory 1221 . 
compared to some estimates based on top-down models [7| (Figure taken from I22I L 


Table 1. Numerical examples for the relation rim ^ F’™”' for the minimum fraction F™’" 
that could be excluded with rim events (or: the minimum number of events ^ Fj 
required to exclude a fraction F^) for a confidence level a = 95%. 


6e^39.3% 10e^25.9% 30e^9.5% 100^^3.0% 300^^1.0% lOOOe^O.3% 


fraction that could be excluded for a given event number n„i (or in turn: the 
minimum event number n™'" required to possibly exclude a certain fraction F-y) is 
given by 

F™ = 1-(1-«)!/"», and = (1) 

with a being the confidence level of rejection. This theoretical limit is reached only 
if for each individual shower, the observations allowed us to completely rule out 
photons as primary particle. Some numerical examples are listed in Tab. Q 

5 Photonuclear cross-section 

The experimental results for UHE photons are based on comparisons to photon 
shower simulations. For the simulation, an extrapolation of the photonuclear cross- 
section is needed. This is in particular important for the production of secondary 
muons in primary photon showers. For hadron primaries, most muons are produced 
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Fig. 6. Data and PDG extrapolation [ssiEni of the photonuclear cross-section 

a^p. Also shown are two parametrizations with larger cross-sections at nltra-high energy, 
denoted 1811 and a'^^° 1271 (see text), and another recent fit 1801 . 


in the decay of pions and kaons generated in collisions of shower hadrons with the 
target air nuclei. In photon-induced showers, however, photoproduction is the main 
process that transfers energy from the electromagnetic to the hadronic channel. 
Muons are then produced in a second step in the hadronic sub-showers that were 
initiated by a photonuclear interaction. Changing the photonuclear cross-section 
influences the rate of transferring energy to hadrons and, thus, the production of 
secondary muons. 

Also the position of depth of shower maximum is affected: Increasing, for in¬ 
stance, the photonuclear cross-section makes primary photon showers more similar 
to hadron-induced cascades. Qualitatively, in this example, the Xmax values would 
become smaller for UHE photon showers, i.e. values closer to the average Xm&x for 
primary hadrons (see Figure P|l would emerge. 

In the extreme case, an unconverted primary photon can undergo photoproduc¬ 
tion in the first interaction of the shower cascade. This would result in a shower 
that hardly can be distinguished from that of a hadron primary. The probability 
for such a process is small, though. The electromagnetic Bethe-Heitler pair pro¬ 
duction cross-section is about 500 mb whereas the photoproduction cross-section is 
expected to be of the order 10 mb at cslO EeV. However, the LPM effect reduces 
the Bethe-Heitler cross-section at ultra-high energy, and an uncertainty exists for 
extrapolating the photonuclear cross-section. 

In FigureEl data and some published extrapolations for the cross-section a^p are 
shown. Large differences between the extrapolations are evident. With respect to 
photon shower studies, several questions arise: Is there a theoretical maximum for 
the allowed increase of the cross-section (smaller values would make the experimen¬ 
tal limits to UHE photons even more severe)? What is the impact on photon shower 
features when adopting different cross-section extrapolations? Is photon shower dis- 
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crimination still possible even for very rapidly rising cross-sections? And, finally: 
If we were to observe photon showers, could this be used in turn to constrain the 
photonuclear cross-section in an energy range inaccessible at colliders? 

To answer these questions, photon shower simulations were performed for vari¬ 
ous cross-section assumptions with CORSIKA taking preshower formation HH] 
and the LPM effect |ld[l24| into account. We consider here the extrapolation of the 
Particle Data Group, ESI 123, as a baseline assumption. The most rapidly 

rising photoproduction cross-section is obtained if the hard pomeron model of Don- 
nachie and Landshoff EZI is extrapolated to high energy. This model does not ac¬ 
count for unitarity corrections that are expected to become important at very high 
energy. Therefore, it should be considered as a rather extreme extrapolation that 
could represent an upper limit of possible low-energy data extrapolations. In the 
following, we will refer to this cross-section extrapolation as (exotic). 

What is the impact on photon showers? The impact due to switching from 
to is quite significant, increasing the number of ground muons by roughly 
70—80% E2- The effect on A^ax is a rather moderate reduction by ~30 g cm“^ 
for photons that formed a preshower before entering the atmosphere EBIj see Fig¬ 
ure 0 The reduction is found to be much larger (up to 100 g cm“^ or more) for 
unconverted photons. 

Recent theoretical work on a possible maximum cross-section seems to indicate 
that cross-section values larger by a factor 2 or more than are very un¬ 

likely EHl. A recent fit to the low-energy data, EH, stays also close to 
Using these cross-sections, the uncertainty of ultra-high energy shower predictions 
is rather moderate. For example, assuming the extrapolation EU (see Fig.|H) 
results in a change of ~7 g cm“^ in Amax and ~10% in muon number. Thus, it 
seems that a reasonable estimate of the present uncertainty from extrapolating the 
photonuclear cross-section is ~10 g cm“^ and ~15% for A„iax and muon number, 
respectively. 

Is photon shower discrimination possible even for rapidly rising cross-sections? 
Assuming an exotic scenario such as tr®™, the discrimination power for photons 
is reduced, but still photons are on average expected to have fewer muons (factor 
2) and larger Amax (50—100 g cm“^) compared to protons. The separation to 
heavier nuclei is even larger. Thus, if UHE photons exist, we have good chances to 
identify them. In particular, the preshower effect can be used as a tool for photon 
identification: The detection of the characteristic dependence of A^ax and muon 
number, or of observables related to these, on primary direction and energy (both 
the average values of the shower observables and their fluctuations) would be an 
unambiguous signal of primary photons. 

Could the photonuclear cross-section be constrained at ultra-high energy by air 
shower observations? When UHE photon showers are identified, their observed 
shower features could in turn be used to constrain the photonuclear cross-section 
at highest energies. The differences in Amax and muon number when assuming 
^PDG Qj. ^exo determined even with relatively small event statistics. Thus, 

it might be possible to convert the observed shower features into an upper limit 
to the photonuclear cross-section. It may also be useful to look for specific photon 
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Fig. 7. Shower maximum distribution of primary photons compared to the reconstructed 
value of the Fly’s Eye event. The measured depth is shown with the la- and 2a- uncer¬ 
tainty. In addition to the simulations with standard photonuclear cross-section, results 
when assuming the extrapolation a‘^^° (see Fig. EJ are given |28|. 


event classes: For instance, photon events with a very delayed first interaction (due 
to LPM suppression) would be much less frequent in the cr®’^° scenario, as they 
would start a hadronic cascade at some point. The rate of deeply starting photon¬ 
like events would constrain the photonuclear cross-section at these primary energies. 
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